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ABSTRACT 

Swift discovered GRB050128 with the Burst Alert Telescope and promptly pointed its narrow field instruments 
to monitor the afterglow. X-ray observations started 108 s after the trigger time. The early decay of the afterglow 
is relatively flat with a temporal decay modeled with a power law with index ~ —0.3. A steepening occurs at 
later times (~ 1500 s) with a power law index of ~ —1.3. During this transition, the observed X-ray spectrum 
does not change. We interpret this behaviour as either an early jet break or evidence for a transition from the fast 
cooling regime to the slow cooling regime in a wind environment. 

Subject headings: gamma rays: bursts; X-rays: individual (GRB050128) 
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1. INTRODUCTION 

The Swift Gamma-ray Burst Explorer (Gehrels et al. 2004) 
was successfully launched on 2004 November 20. Thanks to its 
fast-pointing capabilities, Swift is performing the first compre- 
hensive observations of the early afterglow phase of Gamma- 
Ray Bursts (GRB). A few GRBs have been followed by Swift 
within 200 s from their trigger time: GRB050117a (193 s), 
GRB050126 (131 s), GRB050128 (108 s), GRB050215b (108 
s), GRB050219a (92 s) and GRB050315 (83 s). 

In this paper we focus on GRB050128. The Burst Alert 
Telescope (BAT; Barthelmy et al. 2005) on board Swift trig- 
gered and located GRB050128 at 04:19:54 UT (Cummings et 
al. 2005). The burst profile is multi-peaked with a Tgo du- 
ration of 13.8 s. The fluence is 4.5 x 10~ 6 erg cm -2 (15- 
350 keV) making it a 'normal' burst with respect to the BATSE 
GRB population. The spectrum of the burst during the Tgo in- 
terval can be described by a power law model with photon index 
T = 1.5±0.1(15-350keV).Thepeakenergyisabove350keV 
making it a classical GRB. 

Swift pointed autonomously to the GRB. We will report 
in the next sections about the XRT and UVOT observations. 
Ground-based follow-up observations started as soon as the 
GCN circular announcing the discovery of the new GRB was 
issued. This happened with some delay, since Swift was in 



the early phases of the mission and each circular was being 
checked manually before being distributed. The first GCN cir- 
cular on GRB050128 was issued by the XRT team (Antonelli 
et al. 2005). The robotic 60-cm REM telescope located in La 
Silla pointed to GRB050128 approximately 3 hr after the burst 
with good seeing conditions (~ 1"). No new sources were dis- 
covered with an upper limit of H > 17 (Covino et al. 2005), 
V > 18.2, R > 18.2 and / > 17.9 (Melandri et al. 2005). A 
further upper limit came from the 2-m Faulkes telescope South 
with R > 20.5 1 1.5 hr after the burst (Monfardini et al. 2005). 
GRB050128 has also been observed in the radio band at 8.4 
GHz yielding an upper limit of 100 /iJy ~ 11 d after the burst 
(Frail & Soderberg 2005). 

In the following we focus on the observations by the X-Ray 
Telescope (XRT, Burrows et al. 2005a) on board Swift. In sec- 
tion 2 we describe the data analysis. In Section 3 we discuss on 
theoretical implications of these observations and in Section 4 
we draw our conclusions. 

2. XRT AND UVOT OBSERVATIONS 

UVOT (Roming et al. 2005) observations started on Jan 28, 
2005 at 20:03:07 UT, about 15.8 hours after the burst. This de- 
lay happened because UVOT was in safe mode at the trigger 
time. During each orbit UVOT acquired a single shot image in 
each filter U, B and V in imaging mode. UVOT did not detect 
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the afterglow. Upper limits on the first (summed over about 1 
d elapsed time) UVOT observations are V < 18.9 (V < 19.8), 
B < 20.3 (B < 21.1) and U < 19.2 (U < 20.3). 

XRT observed GRB050128 after an automatic slew of the 
Swift satellite. However, XRT was not yet operating in its 
nominal automatic mode changing configuration but rather in a 
manual mode for the purpose of obtaining calibration data. Be- 
fore pointing to GRB050128, the XRT was observing a UVOT 
calibration target in photon counting mode (see Hill et al. 2004 
for a description of XRT observing modes). XRT fully set- 
tled on the BAT position 108 s after the trigger. XRT observed 
GRB050128 for 17 orbits following the first pointing, accumu- 
lating a nominal exposure time of 17303 s (distributed over 73 
ks). This low Earth orbit of Swift causes source observations to 
be interrupted each orbit. At the same time, thanks to the fast- 
pointing capability of Swift, several targets may be observed 
per orbit. At this early stage of the mission the analysis of the 
data is not straightforward. We analysed the data running the 
task xrtpipeline within FTOOLS v5.3.1 and cutting out 
temporal intervals when the CCD temperature was higher than 
—50° Celsius (see Burrows et al. 2005b) and when the total 
count rate in the 0.2-10 keV energy band over the entire CCD 
was larger than 85 c s _1 (these counts are mainly soft counts 
and are due either to dark current or to the bright Earth limb 
near the end of each snapshot observation). With these cuts 
we obtain a total exposure time of 13047 s distributed over 1 1 
orbits. 

2.1. Angular analysis 

A fading source is clearly evident in all the XRT orbits. In 
the first two orbits the source is clearly piled-up and to derive 
an unbiased position we rely on the remaining <~ 10 ks ex- 
posure. An image has been extracted in the 0.5-10 keV en- 
ergy band to avoid contamination from low energy photons. 
The source position has been derived with XIMAGE (v4.2.1) 
using the cent roid command: RA(J2000): 14h38ml8s.O, 
Dec(J2000): -34°45'55".9. The main contributors to the po- 
sitional uncertainty are the uncalibrated satellite attitude and 
boresight, resulting in a ~ 6" error radius (90% confidence 
level). 

2.2. Temporal analysis 

In order to properly track the decay of the fading source we 
have to account for the piled-up core in the first two snapshot 
observations. To this aim we extracted photons from an annular 
region (inner and outer radii 4 and 30 pixels, respectively) on 
source. This aperture was then applied to the rest of the obser- 
vations, even when it was not needed. The light curve will have 
an underestimated normalization but it will not be distorted by 
pile-up. A background light curve has also been extracted from 
an annular region (inner and outer radii 80 and 120 pixels, 
respectively) centered on the same position and free of other 
sources and hot pixels. Fig. 1 shows the background-subtracted 
light curve in the 0.2-10 keV energy band. The source is clearly 
fading. The decay light curve is not consistent with a single 
power law ix^. ed = 2.6, with 35 degrees of freedom, dof, null 
hypothesis probability, nhp, 7 x 10~ 7 ), being flat at the begin- 
ning and steepening at later times. We also tried a single power 
law with a different initial time to, considered here as a free 
parameter. We can account in principle for the observed decay 

15 See e.g. the XMM-Newton User's Handbook at http://xmm.vilspa.esa.es/extt 



with to = —780 ± 290 s, however this early time is not phys- 
ically meaningful. A better description of the data is provided 
by two power laws (with indices ct\ and a^) smoothly joined at 
a break time 4 to model the decay. The fit with this model is 
good ( X 2 red = 0.7 with 33 dof, 0.90 nhp) and on = -0.27 
a 2 = — 1.30±g;if and t b = 1472±|°^ s (these errors are 90% 
for one interesting parameter, i.e. A% 2 = 2.71, throughout the 
letter). 

2.3. Spectral analysis 

Given the large variability in the source count rate and the 
knowledge from the temporal analysis of the existence of a 
change in the decay slope, we extracted three spectra from 
our data, one from each of the first two snapshots and one for 
the rest of the observation (see Fig. 1). The first two spectra 
were extracted from the same annular region as for the tempo- 
ral analysis. The last spectrum, since the source is much fainter, 
was extracted from a circular region on source of 30 pixel ra- 
dius. Exposure times are 286, 1653 and 10731 s, respectively. 
Data have been filtered for grades 0^- (according to the XMM- 
Newton nomenclature, i.e. single and double pixel events 15 ), 
since at this stage the response matrix for the standard 0-12 
pixel event is not fully calibrated. Ancillary response files were 
generated with the task xrtmkarf within FTOOLS (v5.3.1), 
accounting for the different extraction regions. Data were re- 
binned to have at least 20 counts per energy bin and allow \ 2 
fitting within XSPEC (vl 1.3.1). 

We fit the data with an absorbed power law model for all 
the observations. We first fixed the absorbing column density 
to the Galactic value of N H = 4.8 x 10 20 cm" 2 . We added 
a 5% systematic uncertainty to all our fits to account for the 
residuals still present in our response matrix (given the rela- 
tively low number of counts this does not alter our results sen- 
sibly). This simple fit can account for the observed spectra. 
The fit is good with xled = 11 for 63 dof (°- 21 nh P)- ^ 
power law photon index is T = 1.66^q 07- The 0.2-10 keV 
unabsorbed fluxes of the three observations with mid times of 
258 s, 6156 s and 51587 s are 2.2 x lO" 10 , 1.9 x 10" 11 and 
6.6 x 10~ 13 erg s _1 cm~ 2 , respectively. Given the slope 
change in the light curve we also untie the power law photon 
index of the three observations (even if it is not required by the 
statistics). The three photon indices are 1.59±0.08, 1.79±0.11 
and 1.59 ± 0.23, respectively. These values are consistent with 
each other with a small deviation in the second snapshot. Al- 
though the fit is consistent with the Galactic column density, 
we let the column density value free to vary. The fit is im- 
proved with a xl e d — 0-9 (62 dof, 0.66 nhp, see Fig. 2) and the 
improvement is significant according to the F-test (probability 
4 x 10~ 4 , even if we improved an already statistically good 
fit). In Fig. 3, we show the contour plot of the column density 
vs. power law photon index. The fit obtained with the column 
density fixed to the Galactic value is outside the 3<r bound- 
ary. The absorbing column density is (1.0 ± 0.2) x 10 21 cm~ 2 
and the power law photon index T = 1.88 ± 0.12. Unab- 
sorbed fluxes (0.2-10 keV) are 2.4 x 10~ 10 , 2.0 x 10~ n and 
7.0 x 10 -13 erg s" 1 cm~ 2 , respectively. Also in this case, 
leaving free the photon index to vary within the observations 
the second one is characterized by a slightly steeper index. 

3. PHYSICAL INTERPRETATIONS 
mm_user_support/documentation/uhb/index.html 
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The major result on the GRB050128 afterglow concerns the 
monitoring in the X-ray band of its early temporal decay. This 
decay cannot be described by a simple power law but can be 
accounted for by a slowly varying double power law decay. 
During this transition there are no apparent marked spectral 
changes. The most straightforward interpretation is that the 
temporal break reveals a jet, i.e. corresponding to the epoch 
when the relativistic beaming angle (I/7) becomes larger than 
the physical opening angle (6j) of the jet during the fireball 
deceleration (e.g. Rhoads 1999). In the slow cooling regime, 
for a uniform density circumstellar medium, the temporal de- 
cay changes from i 3 ^-?)/ 4 to - t~ p (e.g. Rhoads 1999), 
which is well consistent with the observed temporal decay in- 
dices when p ~ 1.3 is adopted. In such a case no spectral 
change is expected. However, the expected spectral photon in- 
dex should be — (p + l)/2 ~ —1.15, too small to be com- 
pared with the observed value. In order to make the jet model 
work, one needs to assume v c < vx < v m before the jet 
break, and vx > max(i/ c ,f m ) after the jet break (here v m 
and v c are the typical synchrotron frequency and the cooling 
frequency, respectively). In such a case, p ~ 1.3 gives a con- 
sistent interpretation of both spectral and temporal indices in 
all three epochs, regardless of whether the medium is an in- 
terstellar medium (ISM) or a wind from a massive companion. 
This model requires a little bit of coincidence in that the syn- 
chrotron frequency happens to cross the X-ray band during the 
jet break. However, considering the rapid decline with time of 
v m this is not a very unlikely possibility. Another caveat is that 
a flat electron spectrum p ~ 1.3 is abnormal in late afterglow 
fits (e.g. Panaitescu & Kumar 2001). However, since we are 
observing a previously unexplored early epoch, a small p re- 
quired for the jet model to work cannot be ruled out. Possible 
ways to generate a flat electron spectrum have been suggested 
earlier (e.g. Bykov & Meszaros 1996). If this is indeed a jet 
break, this would be the earliest jet break detection so far. Us- 
ing the standard definition of jet break time (i.e. 9j = l/j(tj)) 
one can derive 8 = 1.8°(i j /2000 s) 3 / 8 {E^/n)- 1 ' 8 [(1 + 
z)/2}~ 3 / 8 for a constant density interstellar medium and 9 = 
3.7°(t,/2000 s) 1 / 4 (£ 52 /A0- 1/4 [(1 + z)/2]- 1 / i for a wind 
model. Here n is the density of the ISM, E52 is the isotropic- 
equivalent burst energy in units of 10 52 erg and A = M/(4ttv) 
is the wind parameter, with M being the mass loss rate, v be- 
ing the wind velocity, and A* = A/(5 x 10 11 g cm -1 ). These 
jets are not extremely narrow (e.g. Covino et al. 2003) but 
are narrower than the typical jets identified in the previous late 
afterglow observations (e.g. Table 2 of Bloom et al. 2003). Ac- 
cording to the GRB standard energy argument (e.g. Frail et al. 
2001 ; Panaitescu & Kumar 2001), such a narrow jet should cor- 
respond to large isotropic gamma-ray energy. Since this burst 
was not particularly bright, it might lie in the low energy tail 
of GRB -energy distribution, thus being another outlier for the 
standard energy relation. 

Besides the jet interpretation, one could search for other pos- 
sible solutions by considering the temporal and spectral rela- 
tions in various afterglow models (e.g. Meszaros, Rees & Wi- 
jers 1998; Sari, Piran & Narayan 1998; Chevalier & Li 2000; 
Zhang & Meszaros 2004). The most straightforward model 
is within the framework of the standard isotropic wind model 
(Chevalier & Li 2000). The first cluster of the data corre- 
sponds to the v c < vx < v m regime, in which the temporal 
index —1/4 and the photon spectral index —3/2 are expected. 
The second and the third clusters of the data correspond to 



the regime of v m < vx < v c , in which the temporal index 
— (2>p — 1 ) /4 and the photon spectral index — 1 — (p — 1 ) /2 are 
expected. All these are consistent with the data for p ~ 2.1. 
In this interpretation, one needs to assume that both v m and 
v c cross the X-ray band during the gap between the first two 
clusters of data and that the frequencies switch the order. This 
could be achieved with a small wind parameter (e.g. A* in the 
range of 0.01 — 0.001). One caveat is that in the wind model 
v c oc t 1 / 2 , so that the time interval of the gap is not long enough 
for v c to completely cross the band. Nonetheless, the spec- 
tral slope in the second cluster is slightly steeper than the other 
two, which might be still consistent with the data if one intro- 
duces an evolving cooling break near the high energy edge of 
the band during the epoch of the second snapshot. Introduc- 
ing a slightly steeper density profile (larger than r~ 2 for the 
wind case) could further alleviate the problem. Furthermore, if 
the electron equipartition factor €3 drops during the temporal 
gap, this would speed up the v c crossing time scale, making the 
model more consistent with the data. 

Finally, we note also that a similar behaviour has been ob- 
served in optical light curve GRB021004 even if with a slightly 
longer break time (~ 0.1 d, Fox et al. 2003). This transition 
has been interpreted by Li & Chevalier (2003) as a fast to slow 
cooling transition. 

4. CONCLUSION 

Swift is exploring for the first time the early stages of GRB 
afterglow decays. We detect a clear early temporal break in 
the X-ray afterglow of GRB050128, with the spectral indices 
not changing around the break. The data could be argued to 
be consistent with either a jet model or a wind model. The jet 
model requires a flat electron spectrum and an assumed spec- 
tral domain change within the temporal gap between the first 
two snap shot observations. If this is true, we may have de- 
tected the earliest jet break so far. The wind model requires a 
(relatively) low ambient density and possibly an evolution of 
the 6b parameter. We note that in this last case, passages from 
fast to slow cooling regimes might have remained hidden in the 
great majority of GRB afterglows if they are characterized by 
such a tenuous environment, due to the late times at which they 
were observed. 

The early detection of the X-ray afterglow, coupled with 
the initial flat decay, allows us to estimate its fluence T = 
J F(t)dt. In fact, given the value of a.\ and a<i, T is not very 
sensitive to the start time t of the afterglow and is dominated 
by the flux at the break time 4. The fluence between 108 s and 
73000 s is (7 ± 2) x 10~ 7 erg cm" 2 (15-350 keV), while ex- 
trapolating from t = to +00 it is T = 9+ 3 x 10~ 7 erg cm~ 2 . 
These values amount to 15% and 20% of the prompt fluence 
in the same energy band. Since prompt and afterglow spectra 
are similar we might expect a relatively small difference in the 
bolometric correction. This is the first determination of the ratio 
between GRB proper and early afterglow energetics. 

This work is supported at OAB by funding from ASI on grant 
number I/R/039/04, at Penn State by NASA contract NAS5- 
00136 and at the University of Leicester by PPARC on grant 
numbers PPA/G/S/00524 and PPA/Z/S/2003/00507. We grate- 
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Time from burst (s) 

FIG. 1. — XRT 0.2-10 keV light curve extracted from an annular region centered on GRB050128. The continuous line represents the fit with two power laws 
smoothly joined. The conversion factor to translate the count rate into a 0.2-10 keV unabsorbed flux is 8.7 X 10~ n erg cm~ 2 counts -1 (for a freely varying 
absorbed power law model, see text). 



6 




0.5 1 2 5 10 

Energy (keV) 



FIG. 2. — XRT 0.2-10 keV energy spectrum of GRB050128. In the upper panel are plotted the spectra of the three snapshot observations described in the text (1: 
squares, 2: circles, other: stars) fit with an (freely) absorbed power law model. In the lower panel there are the residuals from the same power law fit to all the data. 
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FIG. 3. — Contour plot of the column density versus the power law photon index for the X-ray spectrum of GRB050128. Contours refer to 1, 2 and 3cr confidence 
level. At the left of the contour plot the Galactic column density interval is drawn centered on the value of 4.8 X 10 20 cm~ 2 and with a 15% uncertainty. 



